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ABSTRACT 

We use the hot spot echpse times of the newly discovered deeply-echpsing dwarf nova 
lY UMa to trace out the shape of its disc during the late superhump era. We find an 
eccentric disc. We show that the brightness of the stream-disc impact region varies 
as expected with | AT^p, where AV is the relative velocity of the stream with respect 
to the velocity of the disc at the impact point. We conclude that the hot spot is the 
source of late superhump light. 

Key words: accretion discs - binaries: close - binaries: eclipsing - stars: individual: 
lY UMa - stars: cataclysmic variables 



1 INTRODUCTION 

lY UMa was identified as a dwarf nova type cataclysmic 
variable (CV) in January 2000 (Uemura et al. 2000) when 
it went into superoutburst showing strong superhumps and 
deep eclipses in its lightcurve. 

Superhumps are luminosity variations with period, P^h, 
a few percent longer than the orbital period, thought to arise 
from the interaction of the donor star orbit with a slowly 
progradely precessing non-axisymmetric accretion disc. The 
eccentricity of the disc arises because a 3:1 resonance occurs 
between the donor star orbit and motion of matter in the 
outer disc. This can only occur in systems with a sufficiently 
low mass ratio (q = .'1°"°'' ) that the 3:1 resonance radius 
is within the tidal radius at which the disc is truncated by 
tidal forces (Paczynski 1977). The superhump period is then 
given by 
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where Pproc is the disc precession period and Po^b is the 
orbital period. 

A detailed study of the superoutburst of lY UMa was 
carried out by Patterson et al. (2000) (hereafter P2000). lY 
UMa has orbital period Porb = 1.77 hours. The system pa- 
rameters were estimated in P2000 as M^d = O.86±O.11M0, 
Mio-nor = 0.12 ± O.O2M0 and orbital inclination i = 86?8 ± 
1?5. 



2 OBSERVATIONS 

Here we analyse V-light data covering the late superhump 
era, HJD 2451572-79, from the more extensive data pre- 
sented in P2000. These data comprise 22 eclipses with typi- 



Exciinple orbital lightcurve 




Figure 1. An example orbital lightcurve of lY UMa during the 
late superhump era which clearly shows the white dwarf and hot 
spot eclipse ingresses and egresses. This orbit has mid-eclipse at 
HJD 2451578.69. 



cal time resolution 20 to 30 seconds. Figure fy shows a repre- 
sentative orbital lightcurve, where a strong orbital hump 
is visible as the stream disc impact region (bright spot) 
comes into view before eclipse. The white dwarf and hot 
spot eclipses are easily identified. 

The entire dataset being analysed is shown in Figure 
bk, and is discussed in Section BI 
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Eclipse egreya phase on ballistic streani tr'ajectory 



Figure 2. "Shadows" cast by the donor star at the beginning 
and end of hot spot ingress (observer looking from 1 and 2 respec- 
tively) and at start and end of hot spot egress (3 and 4) constrain 
the location of the hot spot in the orbital plane to a diamond 
shape. 



3 TRACING THE ECCENTRIC DISC SHAPE 

3.1 Shadow method 

The eclipse of the prominent hot spot (Figure hh provides an 
opportunity to probe the disc shape. Assuming the surface 
of the donor star is described by its Roche potential surface, 
and if the mass ratio, q, and orbital inclination, i, are known, 
the region of the orbital plane eclipsed by the donor at any 
particular orbital phase can be calculated. Figure H shows 
how these "shadows" cast by the donor star at phases cor- 
responding to the start and end of hot spot eclipse ingress 
and egress constrain the location of the hot spot in the or- 
bital plane to a diamond shaped region. This method has 
been used before, e.g. Wood et al. (1986). The hot spot is 
the point where the accretion stream impacts the edge of 
the disc, so determining the hot spot location determines 
the location of the disc edge. Assuming the disc precesses 
at the beat period between the orbital and superhump peri- 
ods, the geometry shown in Figure was used to trace out 
the outline of the disc edge using those eclipses where both 
ingress and egress of hot spot eclipse could be measured. 



3.2 Stream trajectory method 

For some eclipses, only the egress of hot spot eclipse could be 
measured, and so the technique described above could not be 
employed for those eclipses. However there is an alternative 
technique for determining the location of the hot spot which 
can be employed in this situation (Hessman et al. 1992). 
Since the hot spot lies on the accretion stream, and the 
times of hot spot egress depend on the distance of the hot 
spot from the line of centres of the two stars, a ballistic 
trajectory for the accretion stream can be used to infer the 
position of the hot spot along the stream. Hence the location 
of the hot spot and the disc edge can be determined. Figure 
shows the variation of eclipse egress phase along a ballistic 
stream trajectory. 

The uncertainty in the time of each point of contact 
translates mostly to an error in the radial distance of the hot 
spot from the white dwarf for an assumed q. Azimuthal error 
was estimated by considering three slightly different stream 
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Distance along stream from Ll/a 

Figure 3. Orbital phase of eclipse egress for points on a bal- 
listic stream trajectory calculated assuming the shape of the 
donor star to be described by its Roche potential surface. The 
solid curves represent a stream started from the LI point with 
(Vx,Vy)={—Cs,0) (see Section Bj) while the dotted curves are 
streams started with {Vx, Vy)={—Cs , ±Cs). The crosses are the 
measured egress phases at their inferred positions along the 
stream. 



Geometry of hot-spot eclipses 




Figure 4. The geometry of the system used to study the 
hot spot eclipses. This figure shows the situation at orbital 
phase 0.11, shortly after a hot spot eclipse. The disc coordinates 
{^diac, Ydisc) are centred on the white dwarf. The stream impacts 
the disc edge (thick dot-dashed outline) at the hot spot, at angle 
Oori) to the line of centres of the two stars, which itself makes an 
angle 6orb with the line of sight. The X-axis of the disc makes 
an angle 0^^^,, with the line of sight, and so the position of the 
stream-disc impact in the precessing disc frame makes an angle 
f^disc = f^orb + ^orb ~ ^disc with the disc X-axis. P is approxi- 
mately the angle between the stream and disc velocities at the 
point of impact. 
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trajectories with initial velocities at the inner Lagrange 
point (with axes as shown in Figure 0) of (14, 14)= (—Cs, 0), 
{—Cs, Cs) and {—Cs, —Cs)- Ca is the sound speed at the LI 
point estimated using 



lY UMa - Disc orientation at late-superliump phase 



Cs ~ 10 ms 



T 



W^K 



where the photospheric temperature T comes from the stel- 
lar models by Baraffe (Baraffe et al. 1998). The diflFerent 
hot spot positions coming from these three slightly diflFerent 
trajectories provide the small near-azimuthal error bars on 
the egress and ingress locations in Figure H. 

3.3 Combining the methods 

The disc radius determined using the stream trajectory 
method depends on whether the start, middle or end of 
egress is used. We find that for those eclipses where both 
ingress and egress were measured, using the mid-egress times 
to locate the hot spot on the accretion stream gives results 
consistent with the robust "shadow" method whose results 
completely constrain the location of the hot spot in the or- 
bital plane. We therefore use the shadow method to deter- 
mine hot spot location and disc radius for those eclipses 
where both hot spot ingress and egress are measurable, and 
the stream trajectory method with mid-egress time for the 
remaining eclipses. 

We use the white dwarf mid-eclipse ephemeris 
Tmw = HJD 2451570.85376(2) + 0.07390906(7)^ 
and late superhump maximum ephemeris Tmid = 
HJD 2451571.731(3) + 0.07558(6)£' from P2000 to cal- 
culate orbital phase and superhump phase. The beat phase 
between orbital and superhump phase, interpreted as disc 
precession phase <j)prec, is defined as zero when superhump 
maximum and mid-eclipse coincide. The angle 9disc (see Fig- 
ure m is then 2-K<pprec- 



4 RESULTS 

The average disc radius Rdtsc in units of orbital separation 
a is < Rdisc/i >= 0.32 with a spread of 0.08. The disc 
shape is shown in Figure ra and is clearly not centred on the 
white dwarf. The shape is similar to that found in OY Car 
by Hessman et al. (1992). 

A good fit to the disc shape was achieved using an ellipse 
with one focus at the white dwarf, semi-major axis adisc and 
eccentricity e, described by 



^disc 



Ojdisc \^ 
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ao is the angle in the precessing disc frame correspond- 
ing to maximum disc radius. The resulting fit has adiac = 
0.34a ± 0.06a, e = 0.31 ± 0.17 and oq = 117 ± 33° and is 
shown in Figure H. The errors in our timing determinations 
and hence deduced geometry (particularly when using the 
stream trajectory method) are rather pessimistic; our deter- 
minations of adisc, e and ao are robust to adopting instead 
more optimistic error estimates. 

The disc is within the tidal cut-off radius r ~ 0.9-Rl = 
0.50a. The largest radius part reaches the location of the 3:1 
resonance at r ~ 0.46a (Osaki 1996), providing a possible 




Figure 5. The disc shape traced out from the hot spot ingress 
and egress times. The disc orientation shown is that which would 
occur at superhump maximum. Thin regions with large radial ex- 
tent were determined using the stream trajectory method (Sec- 
tion 3.2). The other regions, which are azimuthally extended (due 
to the azimuthal extent of the hot spot region and the relative 
motion of hot spot and disc frame during hot spo t ec lipse) were 
determined using the "shadow" method (Section 3.1). The dot- 



dashed outline is the fitted disc shape and the dotted circle has 
radius rdrc (the circularization radius for mass ratio q = 0.14) 
and is plotted to make the non-axisymmetry of the disc easier to 
see. 



explanation for the maintenance of disc eccentricity at this 
stage after the superout burst. The smallest radius coincides 
with the circularization radius. 



5 DISCUSSION 

The disc shape we find is similar to that found for OY Car 
by Hessman et al. (1992), except that we find the minimum 
radius at disc azimuth adiac ~ —63° while in OY Car it is 
around Udiac = 0. This means the smallest radius of the disc 
is lined up with the donor star at superhump phase -0.18 in 
lY UMa while this occurs close to superhump phase in OY 
Car. This has implications for the viability of the hot spot 
as the source of the superhump light during the late super- 
hump era. Superhumps caused by the variation in hot spot 
brightness as the stream impacts the disc at varying depth in 
the white dwarf potential well were first suggested by Vogt 
(1981). This theory has generally been rejected as an ex- 
planation for the common superhumps occurring during the 
decline from superoutburst maximum in SU UMa stars, but 
the possible link between this model and late superhumps 
(which appear during very late decline of the superoutburst) 
has been suggested several times, e.g. by Whitehurst (1988) 
and by Rolfe, Haswell & Patterson (2000) (hereafter RHP) 
after a study of persistent superhumps in the novalike V348 
Pup. See RHP for more discussion. 
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Measured disc radius and fit 




Figure 6. The disc radius as a function of adisc (see defini- 
tion in Figure Kj) . Points marked witli tria ngles were determined 
using the "shadow" method (Section 3.1), while points marked 



with sta rs w ere determined using the stream trajectory method 
(Section Eih. 
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Figure 7. (a) The histogram shows the peak bright spot mod- 
ulation, obtained by subtracting the normalized average orbital 
curve from the normalized data and then extracting the data 
for the orbital phase range 0.85—0.95 and folding on superhump 
phase, (b) shows the measured hot spot eclipse egress intensity 
A//i3e (crosses). The dashed curve in (a) and (b) is the kinetic 
energy of the stream at the hot spot, while the dotted curve is 
the kinetic energy of relative motion of the stream and disc at the 
hot spot. The solid curves are the predicted hot spot peak and 
egress intensities combining our simple 3D model for the hot spot 
structure with the relative KE model for hot spot brightness. See 
Section H The latter three curves were predicted from the fitted 
disc shape and scaled in flux to have the same maximum as the 
measured data, (c) shows the angle /3 between the stream and 
disc velocities at the impact point. 



In Figure Wa we show a measure of the orbital hump 
peak modulation while in Figure Mb we show the hot spot 
flux uncovered on egress of the hot spot eclipse (A/hse from 
P2000). The variation in hot spot brightness, Fhotspot, pre- 
dicted by our fitted disc shape is shown in Figures Ma & b, 
calculated using three models described below. 

(i) The dashed curve in Figure [7k & b is calculated using 
Fhotspot 'X 1/r where r is the radius of the disc at the hot 
spot. This assumes the energy released at the hot spot varies 
as its depth in the white dwarf gravitational potential. 

(ii) The maximum energy which could be released at im- 
pact is the kinetic energy of the relative motion of the stream 
and disc at the hot spot. If the stream velocity vector at the 
point of impact is Vstream and the disc velocity at this point 
is Vdisc then the kinetic energy per unit mass of the relative 

velocity is KE = ^{Vdisc ~ Vstream) -{Vdisc — Vstrcam) ■ We 

approximate 



Vdisc I — \j GM^d ( 



at the impact point with direction parallel to the disc edge 
in the white dwarf frame. This is the correct velocity for an 
elliptical orbit around a point mass. This model is plotted 
as the dotted curve in Figures Ma & b. 

(iii) The solid curves in Figures fm & b were produced 
by assuming the total hot spot fiux behaves as described in 
the previous model, and considering a simple model for the 
3 dimensional structure of the hot spot, shown in Figure pi 
The hot spot has an elliptical cross section in the r — z plane 
with axis of size Vspot in the radial direction and hspot in the 
vertical direction (see inset in Figure H) . rspot and hspot and 
the spot brightness decrease downstream from the initial 
impact as e ' =p°' . Upstream of the impact the hot spot 
surface is rounded off with a hemisphere of uniform bright- 
ness equal to that at the initial impact. The angular extent 
of the hot spot region is set by AOspot ~ A9oadisc/ro, where 
ro is the disc radius at the impact point. This keeps the arc 
length of the hot spot region roughly constant as the eccen- 
tric disc precesses. The disc thickness is 2Hdisc = 0.02a and 
any region of the hot spot surface within the disc is con- 
sidered to be obscured completely. The total flux coming 
from an area element on the hot spot surface is foreshort- 
ened according to its area projected in the direction of the 
observer. This simple parametrized model for the hot spot 
structure enables us to model the full hot spot lightcurve, 
taking into account the visibility of the hot spot at any or- 
bital phase and disc precession phase. The model lightcurve 
is shown in Figure Hb. We use A^o ~ 36° which is the av- 
erage azimuthal extent of the hot spot regions det ermined 
from the eclipse timings as described in Section 



3.1. Noting 



that we see no eclipse of the white dwarf by the hot spot 
region, we determine the maximum possible value of hspot 
for i = 86?8 ± 1?5 to be 0.013a ± 0.006a. Therefore we use 
TspotiP = 0) = hspot{0 = 0) = 0.013a at the point of impact. 

The phasing of the hot spot brightness predicted by 
model (i) does not agree with either of the two measures of 
hot spot variation in Figures |ji & b, with the peak in pre- 
dicted brightness being about 0.25 earlier than the measured 
peak. The fractional amplitude in the predicted variation is 
about half that of the measured variation: in Figures Ma & b 
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Figure 8. The assumed 3D geometry of the hot spot seen looking 
down onto the orbital plane in the —z direction. The hot spot is 
split up into "rings" in the rz-plane which are ellipses centred on 
the elliptical disc edge with axis in the r direction of Vspot and 
vertical axis hspot- Both Vspot, hspot and the hot spot brightness 



decrease as e 



-fl^/A9; 



p°* as the stream merges with the disc, and 



the stream impact is rounded off with a hemisphere for 6 < 0. 



the predicted curve has simply been scaled to have the same 
maximum value as the observed hot spot curves. 

In Figure |7|c we also show the variation of the angle /3 
between the stream and disc velocities at the point of im- 
pact (approximately equal to /3 in Figure y). The hot spot 
emission should arise primarily from dissipation of kinetic 
energy as the infalling stream impacts the disc edge and 
merges with the Keplerian motion of the disc. The available 
kinetic energy increases as the angle /3 between the veloc- 
ity vectors of the stream and disc flows increases. Figure Wc. 
clearly shows /3 to be greatest around the observed maxi- 
mum of the hot spot intensity, which coincides with super- 
hump maximum. This effect is taken into account in model 
(ii) where the maximum kinetic energy available at the im- 
pact point is calculated and depends on the directions and 
magnitudes of the stream and disc velocities. The shape, 
phasing and fractional amplitude of the curve for model (ii) 
shows much better agreement with the hot spot brightness 
measured in Figure 7a. The predicted curve also agrees well 
in phase and fractional amplitude with the measure of hot 
spot flux in Figure Mb, although the hot spot flux between 
superhump phase ~0.6 and 0.9 is much lower than predicted. 

Our final model, (iii), which combines the relative ki- 
netic energy model with a treatment of the hot spot visibil- 
ity, does not differ from the previous model in its prediction 
for the peak flux of the orbital hump (Figure Ra). This is to 
be expected since at the peak of the orbital hump we should 
be seeing the whole length of the bright spot clearly, so the 
visibility at this orbital phase should not be very sensitive 
to small changes in orientation through the disc precession 
cycle. The hot spot flux measured using the hot spot eclipse 



egress is different from the predicted variation in intrinsic 
hot spot flux. This is because at the orbital phase of hot 
spot egress, we are seeing the impact region roughly end- 
on, and at this orientation the fraction of the hot spot flux 
reaching the observer is sensitive to the exact orientation 
of the hot spot. The result is that the predicted variation 
(solid curve in Figure up) rises to maximum later than the 
intrinsic variation, but decreases at about the same time as 
the intrinsic variation. This produces a more flat-bottomed 
curve which agrees better in shape with the observed egress 
flux than models (i) and (ii). The flux around egress flux 
minimum is too low in our models, but this is not surpris- 
ing; our model for the structure of the impact region is a 
very simple one. It is also possible that the hot spot emis- 
sion is strongly anisotropic due to the complicated structure 
of shocks and contact discontinuities we expect to be present 
at the impact. 

Figure Ha shows the evolution of the lY UMa lightcurve 
throughout the late superhump era. The important features 
to note are how strongly the amplitude of the orbital hump 
varies, and how the late superhump is very weak, in many 
cases undetectable, when it does not occur near the orbital 
hump. Strong variation in hump amplitude was also seen 
towards the end of a more recent superoutburst of lY UMa 
(Uemura 2000). The combined lightcurve of the late super- 
hump and the orbital hump is not simply a sum of the two, 
it is some non-linear combination, with the combined hump 
being extremely strong when orbital hump and late super- 
hump coincide, while both orbital hump and late superhump 
are weak when they are not coincident. This is exactly what 
we expect if the stream-disc impact is the source of the late 
superhump. The hot spot is modulated at the late super- 
hump period, as the impact geometry of the stream and 
disc varies. Consequently, the orbital hump, which is just the 
emission from the hot spot coming into view, is modulated 
on the disc precession period, being strongest at late su- 
perhump maximum. The late superhump is very weak away 
from the orbital hump since away from the orbital hump the 
bright spot is on the far side of the disc and hence has a low 
visibility. Our model hot spot lightcurves produced using 
model (iii) are shown in Figure Hb. Ignoring the stochas- 
tic flickering in the observations (which is strongest in the 
second panel, HJD 72.7), and the deep white dwarf eclipse 
in the observations not included in the model, our simple 
model does a convincing job of reproducing the observed 
bright spot behaviour, although the peak of the orbital hump 
when it occurs close to late superhump maximum is sharper 
in the observations than in our model. 

The intrinsic bright spot flux in the lightcurve in the 
third panel (HJD 73.7) appears to be lower than the ob- 
served flux, which make is difficult to see how the shape 
of this curve agrees with observations. Figure He shows the 
model and observations for this curve, with the model scaled 
by a factor of 5. The good agreement between the shape of 
the model and the observations is clear, particularly the dou- 
ble humped nature of each orbit curve and the flattish top to 
the orbital hump before some of the eclipses. More detailed 
modeling is necessary to explain the low intrinsic brightness 
of the hot spot in this region of the disc precession cycle. The 
relative height of the hot spot peak above the disc needed a 
bit of fine-tuning to achieve this match between the model 
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(c) Observations and re — scaled model for HJD 73.7 
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Figure 9. (a) shows all 7 lightcurves from the late superhump dataset analysed in this paper. The date of the start of each lightcurve 
is shown in the top left of each plot as HJD-2451500. Arrows indicate times of late superhump maximum according to the ephemeris 
in P2000. Note the significant variation in the amplitude of the orbital hump over the 7 days, (b) Shows our model of the bright spot 
lightcurve for the same time series as (a). The model assumes the intrinsic hot spot brightness to vary as the kinetic energy of relative 
motion of stream and disc at impact, while a simple model for the 3D structure of the hot spot enables us to predict its visibility as a 
function of orbital phase and precession phase. See Section H. (c) shows the model data from from the 3rd panel of (b) scaled in flux 
by a factor of 5 and shifted up by 1.4 flux units, along with the corresponding observed data. This shows how the shape of the model 
lightcurve is morphologically close to that of the observations, although the intrinsic hot spot brightness predicted by our model is too 
low at the disc precession phase. Note that our model does not include the deep white dwarf eclipse seen in the observations. 



and observed lightcurves, and rspot{0 = 0) had to be less 
than about twice hspot{d = 0). 

It is difficult to avoid the conclusion that the late- 
superhump light in lY UMa is predominantly coming from 
the hot spot, whose brightness is modulated at the super- 
hump period as a result of the varying impact geometry 
of the stream and the disc. The disc shape and orientation 



found in OY Car (Hessman et al. 1992) suggests that the 
late superhumps in that system may also originate from the 
hot spot. An analysis of the OY Car observations like that 
in our Figure M would be valuable. Murray (1996) performed 
SPH simulations which support the hot spot as the source 
of late superhump light, but radiative processes were not 
modelled. More detailed hydrodynamic simulations concen- 
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trating on the stream-disc impact region have been carried 
out (e.g. Armitage & Livio 1998), but there heis been no 
study of the hot spot in such detail for non-circular discs. 
Such work would be valuable in the light of these observa- 
tions. 

We haved assumed the shape of the disc to be fixed 
and constant in the precessing disc frame, while our obser- 
vations of V348 Pup (RHP) and SPH simulations (Haswell, 
King, Murray & Charles 2000) suggest the disc shape and 
size changes during the superhump cycle. The disc shape we 
find represents the disc radius seen by the accretion stream 
as a function of superhump phase, and any effect of chang- 
ing shape should have little effect on the angle /3 at which 
the stream impacts the disc since the shape is traced out 
between radius measurements which are closely spaced in 
superhump phase. 

The transition from the common superhumps which ap- 
pear during the decline from maximum of the superoutburst 
to the late superhumps during the very late decline of the 
superoutburst is easy to explain using a minor modification 
to the tidal-thermal- instability model of Osaki (1989) pro- 
posed by Hellier (2000). Hellier suggests that the tidal and 
thermal instabilities are uncoupled: the disc can drop out 
of the thermal high state at the end of the superoutburst 
while remaining tidally eccentric. At this point, the disc vis- 
cosity dramatically reduces and so the viscous dissipation 
due to tidal stressing which causes normal superhumps will 
be similarly reduced, allowing the late superhumps which 
were drowned out during the superoutburst to be seen. 



6 SUMMARY 
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The timing of hot spot eclipses during the late super- 
hump era reveals the disc shape, which proves to be non- 
axisymmetric. 

The inferred shape and orientation of the disc and a sim- 
ple treatment of the kinetic energy available at the stream- 
disc impact enable us to calculate the hot spot brightness 
throughout the superhump cycle. The resulting variation in 
predicted hot spot brightness matches the observed varia- 
tion well, and is in very close agreement with the superhump 
modulation. Combining these results with a simple model for 
the 3D structure of the hot spot enables our model to re- 
produce the observed lightcurves during the late superhump 
era. 

We conclude that the stream-disc impact region is the 
source of late superhump light in lY UMa. 
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